Breakdown of universality in transitions to spatio-temporal chaos. 
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In this Letter we show that the transition from laminar to active behavior in extended chaotic 
systems can vary from a continuous transition in the universality class of Directed Percolation with 
infinitely many absorbing states to what appears as a first order transition. The latter occurs when 
finite lifetime non-chaotic structures, called "solitons", dominate the dynamics. We illustrate this 
scenario in an extension of the deterministic Chate-Manneville coupled map lattice model and in a 
soliton including variant of the stochastic Domany-Kinzel cellular automaton. 
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The nature of transitions in extended deterministic 
dynamical systems is not very well understood. Since 
few analytical methods are available for such systems, 
it is tempting to map them to stochastic models [jp. 
In a study of the deterministic "damped Kuramoto- 
Sivashinsky equation" Chate and Manneville || in- 
troduced the notion of a universal transition to turbu- 
lence via "spatio-temporal intermittency" jl]|5] . Pomeau 
gave general arguments Q for this transition to be in 
the universality class of Directed Percolation (DP) ^,|). 
These arguments rested on earlier work by Grassberger 
H and Janssen [ fTo| , who conjectured that any stochastic 
process with a unique absorbing state should be in the 
class of DP. To check whether deterministic models could 
be in the DP class, Chate and Manneville introduced a 
very simple coupled map lattice (CML), with the local 
map either performing "laminar" or chaotic motion. 

Surprisingly, the critical exponents of this CML were 
not those of DP [^l) ; in fact they appeared to vary contin- 
uously along the critical line. Grassberger and Schreiber 
jjJI pointed out that the presence of long lived "soli- 
tons" (local excitations that propagate with unit ve- 
locity through the lattice, see Fig 0b-c and may 
lead to long crossover times, and conjectured that the 
true asymptotic behavior of the Chate-Manneville model 
would be in the DP universality class. 

Our central claim is that even solitons with a finite life- 
time can completely change the nature of the transition. 
We base this claim on studies of extensions of (i) the 
Chate-Manneville model || and (ii) the Domany-Kinzel 
stochastic cellular automaton E| , where both extensions 
facilitate the tuning of solitonic properties. In the regime 
of short soliton lifetimes the bulk exponents of our CML 
are consistent with DP; the spreading exponents indicate 
that the CML falls in the universality class of DP with 
infinitely many absorbing states. For both our CML and 
our stochastic model the transition between active and 
inactive states looses its continuous nature and appears 
to become first order for large soliton lifetimes. We shall 
argue that this discontinuity emerges because pairs of 
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FIG. 1. Spacetime plots of our coupled map lattice -(||) 
where inactive (chaotic) sites are white (black) for r~3 near 
criticality, illustrating the increasing role of solitons as func- 
tion of our parameter b: (a) b = 0.2, e = 0.373. (b) b — 0, 
(Chate-Manneville) e = 0.357. (c) b=-0.1,e=0.352. 

solitons can generate new activity upon collisions (see 
Fig. |l|c and Fig. [i]a) while individual solitons do not lead 
to new turbulent activity. We study this phenomenon in 
detail in our stochastic cellular automaton and show that 
the non-universality seen in earlier work is likely due to 
the proximity of this (quasi) first order transition: our 
stochastic model shows transient behaviour which can 
be fitted quite convincingly to non-universal power laws. 
Here we focus on the broad picture; a detailed study will 
be presented elsewhere ]l4[ ]. 

Coupled Map Lattice - The model introduced by 
Chate and Manneville IO] is a ID coupled map lattice 



Ui(n + 1) = f(ui(n)) + -A f Ui(n) 



(1) 



where AfUi(n) = 
When u < 1, / 
f(u) = r(| - \u - 



f(ui-i(n)) - 2f(ui(n)) + f(u i+1 (n)). 
is a standard tent map of the form 
- ||) and u displays chaotic behavior, 
while when u > 1 , / is simply the identity and u displays 
laminar behavior. States in which all sites are laminar re- 
main so, but chaotic sites can "infect" their neighbors due 
to the spatial coupling. The effectiveness of this spread- 
ing of the chaos depends on r and e. Taking the density 
of chaotic sites or "activity" m as an order parameter, 
transitions from a "laminar" state (in which m decays 
to zero) to a "turbulent" state (where m reaches a finite 
value in an infinite system) can be defined. To relate 



1 




FIG. 2. Activity in the model (|2j)-(^P for a systemsize of 
2048 and ensemble of 128 runs, (a) Activity at t = 1000 (white 
= inactive) . (b) Activity as function of 8e (distance to critical 
point) at time 2.10 5 for b=— 0.1 (squares), fo = (open circles) 
and 6 = 0.2 (closed circles), (c) Illustration of the sharpening 
of the b= — 0.1 transition for increasing times 5.10 3 (+), 5.10 4 
(*) and 5.10 (X); note the strongly magnified scale. 

such CML's to physical systems of weakly coupled 
chaotic elements, we must interpret the map / as a re- 
turn map on a Poincare section; each local map should 
therefore be at least two dimensional and invertible [15| . 
We thus introduce the CML 

Ui(n + 1) = f(ui(n)) + |A/Ui(n) + «j(n) (2) 
v t (n + I) = b(ui(n + 1) - Ui(n)) (3) 

where / is the same as before and the new parameter 
b is the Jacobian of the local map at each site; this 
map is invertible for any non-zero b and becomes in- 
creasingly two-dimensional with |6|. This design is anal- 
ogous to the construction of the Henon map |L6| from 
the logistic map, except that on the right hand side 
b(ui(n + 1) — Ui{n)) appears instead of bui{n). This en- 
sures that the absorbing state fixed points Ui(n) = u* 
of the old CML ([[]) are mapped to the laminar fixed 
point (ui(n) , Vi(n)) = (it*,0). The model is updated 
synchronously and typical dynamical states are shown 
in Fig. [IJ in what follows we will fix the parameter r = 3. 

To get a feeling for the location and nature of the tran- 
sition we show in Fig. ||a the activity after 1000 time- 
steps as a function of b and e; the active and inactive 
phases can be clearly distinguished. In Fig. ||b we show 
that when b is varied, there is a qualitative difference in 
the activity as a function of e. The behavior for b = 
is consistent with a continuous transition, whereas for 
b =—0.1 a marked steepening occurs consistent with the 
emergence of a discontinuity; further support for this 
comes from a study of the correlation function which 
shows a finite correlation length at the transition |fl4) . 
In Fig. ^c we show the sharpening of this discontinuity 
as function of the total integration time. 

CML near continuous transition - We focus now on 
the soliton-poor regime where the transition appears to 
be continuous. We have computed critical exponents us- 
ing finite size scaling pjJTTl]. The critical line in the pa- 
rameter plane (e, r) is located by measuring the "absorp- 
tion time" r(r, e, L), i.e., the average time it takes the 
system, starting from a random initial state, to reach the 
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Table 1. The critical exponents z and 9 = P/vn for our CML 
(§)"(!)■ Note that e c (b) is a multiple- valued function and 
the values for b > —0.1 correspond to the "upper branch" of 
Fig. ^a. The values for DP (last row) are taken from fj"j|| , 

absorbing state. At the critical point e = e c (r), this time 
diverges like r(e c ,L) ~ L z where the usual dynamical 
exponent z = v\\/v\_ has been introduced. The order pa- 
rameter m(e, L, t) is the fraction of chaotic sites in the 
lattice, again averaged over many different initial states. 
The scaling of the order parameter, m ~ (e — e c )^ for 
e — > e+ defines the critical exponent (3. Precisely at the 
critical point the order parameter decays as m(e c , t, L) ~ 
L~ l3 /" ± g(t/L z ). For t < L z this behavior should be in- 
dependent of L and m(e c ,t,L) ~ t~ e with 8 = /3/un. 
Table 1 shows the exponents z and 9 as function of b for 
r = 3. These values where obtained for systemsizes up 
to L = 2048 (128 realizations) H). There are apparently 
regimes for |6| > 0.1 where the exponents are very close 
to their DP values. 

Spreading Properties of our CML - Instead of follow- 
ing the decay of an initially uniformely filled system, one 
can also determine critical exponents from the spreading 
of an initial seed of turbulence in an otherwise laminar 
configuration (T^|. These so-called dynamical exponents 
are defined in terms of the number of chaotic sites N(t), 
the survival probability P(t), the mean-squared deviation 
R 2 (t) from the origin of the turbulent activity and the 
density n(t) of chaotic sites within the spreading patch 
of turbulence. Thus we assume 

N(t) - f P(t) - t- s R 2 {t) - t z ° n(t) - r 6 ' (4) 

For DP one finds ^ = = 6* s and z s = 2/z. In the Chate- 
Manncvillc model at r — 3 it is basically impossible to 
determine these exponents, since the spreading is com- 
pletely dominated by the solitons (Fig. ^b). In our gen- 
eralized CML it turns out to be possible to determine the 
spreading exponents (Fig. ^c) in the weak soliton regime. 
In the CML the absorbing state is non-unique since any 
state in which all it-values are above unity and all v- 
values are not too large is absorbing. We indeed found 
the dynamical exponents to depend on the configuration 
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FIG. 3. Spreading of activity in our CML averaged over 
10 4 realizations for 6= -0.1 (a), 6=-0 (b) and fe = 0.2 (c). 

of the laminar sites; the spreading exponents are non- 
universal. This is consistent with recent studies of sys- 
tems with infinite numbers of absorbing states Jl8| , p0| -p^] 
where it has been revealed that these differ from classi- 
cal DP precisely in the non-universality of the spreading 
exponents. It has also been conjectured, and verified nu- 
merically [ pl| , p2[ , that the dynamical exponents satisfy 
the hyperscaling relation r/ + S + 6 S = dz s /2 where d is 
the spatial dimension. The spreading exponents depend 
on the mean value of the field < u > in the laminar state 
outside of the seed. For the example shown in Fig. |]c, 
the value of < u > which would result from a long run 
with homogeneously random initial conditions is 1.235. 
Even a 2% variation in < u > creates a 100% variation 
of T] whereas A = z s /2 — (77 + 6 + 6 S ) remains less than 
0.1 (for times up to 2000). Similarly, within this range 5 
varies between and 0.5. For other parameter values rj 
may range from 0.8 to 0.1; for more details, see juj. 

First order behavior and Stochastic Model - The cause 
of the discontinuity in the activity in the soliton rich 
regime is illustrated in spacetime plots like Fig. |l|c: while 
individual turbulent patches clearly have a finite life- 
time, new activity is created by collisions of the solitons 
that surround the active patches. To illustrate this sce- 
nario, we extend the Domany-Kinzel cellular automaton 
p3| in the following way: (i) Our model contains two 
species representing the chaotic sites and the solitons. 
(ii) The chaotic sites behave like active sites in usual 
bond-directed DP, except that when they "die" they can 
emit a left or right-moving soliton with probability c. 
(Hi) The solitons travel ballistically and die with a prob- 
ability d. (iv) Individual solitons are seen by the chaotic 
sites as inactive, (v) Upon collision two solitons gener- 
ate a chaotic site; this is the only way by which solitons 
enhance the activity. A typical spacetime plot of our 
model is shown in Fig. |Ja, showing the same qualitative 
behavior as the CML in the soliton rich regime: finite 
size clusters of activity surrounded by clouds of solitons. 

Mean field equations - The rate equations for our 
stochastic model are 



S = cC - S 2 - dS 
j2 „.r<l 



(5) 

C = rC + S 2 - uC' z (6) 
where C and S are the densities of chaotic sites and 



of solitons. Pure DP corresponds to the equation C = 
rC~uC 2 . These mean- field equations display both a 
first and second order transition, and the nature of the 
transition is governed by z=(d/c) 2 u; the change to first 
order behavior occurs at z = l. If z>l (short soliton life 
time) the solitons just renormalize u to u(l — -). When 
z < 1 the behavior for positive r (the active phase) is gov- 
erned by the stable node at S$ ~ a(z" 1 / 2 — 1). When r 
becomes negative this fixed point remains stable and at a 
finite distance away from the origin; simultaneously, the 
origin becomes attractive while a saddle at S = S* « jzj 
appears close to it. Initial conditions close to the origin 
will flow there, but initial conditions above the saddle 
stable manifold will flow to the node. This will go on 
until (S*,C*) and (S£ , C|) mer g e m a saddle-node bifur- 
cation at r = r c {z) < below which the origin becomes 
globally attractive - clearly a first order scenario. 

Activity in first order regime - In Fig. ^b we show 
the evolution of the activity m for our stochastic model 
in the soliton rich regime. There are two important fea- 
tures that can be extracted from this data: (i) For a 
transient period that goes up to time 10 3 , it is possible 
to find values of p such that the decay of m appears to 
be a powerlaw with a non-DP exponent. For the exam- 
ple shown, a reasonable scaling can be obtained over 2 
decades. However, for this to be real asymptotic scaling, 
one should be able to have this scaling extend to arbi- 
trary large times; however the activity curves for sizes 
200, 2000 and 20000 precisely curve down at the same 
time (Fig. ^c); hence there is no hope that increasing the 
systemsize extends the time interval over which apparant 
scaling can be found, (ii) For long times the activity ei- 
ther decays rapidly, or first hits a plateaux. Clearly, for 
increasing time the curve of m as function of p will make 
a sharper and sharper jump, similar to what we found for 
the CML. In the plateaux regime, the qualitative dynam- 
ics is as shown in Fig. |Ja. We have also checked that for 
p = 0.621 the same plateaux value is reached for initial 
activities in the range from 1 to 0.1; for initial activi- 
ties of 0.05 and smaller, there is an initial increase of the 
activity but the plateaux is never reached 0] . 

It remains to be understood whether this soliton- 
assisted turbulent state can persist forever and thus 
whether the observed first order transition is a true ther- 
modynamic transition. One can argue that below the DP 
threshold (i.e., for negative r in (0)) the active state can 
be destroyed by creating a sufficiently large laminar hole, 
so large that the solitons emitted from the edges cannot 
penetrate through it and collide to create turbulence [E5| . 
Thus the motion of the edges of such a large "droplet" 
would only be driven by the diffusion of the chaotic sites, 
which for negative r will cause the laminar region to ex- 
pand. Something like this is observed for states close 
to but below the observed transition, where the activity 
remains almost constant for a long transient time after 
which it rapidly decays. We have not been able to 
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FIG. 4. (a) Soliton dominated dynamics in the stochastic 
model for d — 0.01, c = 0.3 and p = 0.55. (b) Average activity 
m for c = 0.1, d = 0.01 and p = 0.612, . . . , 0.621 for an en- 
semble of 20 systems of size 2 x 10 4 . The dotted curves are 
for comparison and have slopes 0.159 and 0.24 (c) Compar- 
ison of activity for systemsizes 200 (2000 realizations, dots), 
2000 (200 realizations, dashed) and 20000 (20 realizations, 
full curves) for p = 0.613, 0.615, 0.617, 0.619 and 0.621. For 
smaller systemsizes the final fall-off occuring for large p seems 
more gentle, but clearly the early time behavior and the value 
of the plateaux activity are insensitive to finite size effects. 



determine if this transient time actually diverges at the 
apparent first order transition point (which would make 
it a bona fide phase transition). Such an investigation 
is difficult since for finite lattices the true asymptotic is 
always the absorbing state. 

Discussion The overall picture that emerges from our 
models is that the transition to spatio-temporal intermit- 
tency is strongly influenced by coherently traveling "soli- 
tons" , which, even though they have a finite life time, 
change the nature of the transition. Recent work p(| 
showing that the transition in the damped Kuramoto- 
Sivashinsky equation appears first order provides ad- 
ditional support for this scenario. Further support 
comes indirectly from the finding |27j that the Chate- 
Manneville model yields critical behavior consistent with 
DP when it is driven asynchronously, which destroys 
the solitons. We believe that our findings also explain 
the nonuniversality observed in these systems, since our 
stochastic model shows transient behaviour which can be 
fitted convincingly to power laws that vary strongly with 
the parameters |Q . In general our study highlights that 
it is far from trivial to decide which are the important de- 
grees of freedom when mapping deterministic to stochas- 
tic behavior: apparantly innoculous structures may have 
an unexpected strong effect on the coarse grained dy- 
namics and render the "natural" , most simple stochastic 
models inapplicable. 
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